Fluctuations in viscous fingering 
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Our experiments on viscous (Saffman- Taylor) fingering in Hele-Shaw channels reveal finger width fluctuations that were 
not observed in previous experiments, which had lower aspect ratios and higher capillary numbers Ca. These fluctuations 
intermittently narrow the finger from its expected width. The magnitude of these fluctuations is described by a power law, 
which holds for all aspect ratios studied up to the onset of tip instabilities. Further, for large aspect ratios, the mean 
finger width exhibits a maximum as Ca is decreased instead of the predicted monotonic increase. 
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When a less viscous fluid displaces a more viscous fluid 
in a Hele-Shaw channel (a quasi-2D geometry in which 
the width w is much greater than the channel thick- 
ness 6), the interface between the fluids forms a pat- 
tern of growing "fingers" . A single finger forms at low 
flow rates; more complex branched patterns evolve at 
high flow rates. This phenomenon is the prototypical 
example of moving interface problems, such as dendritic 
growth and flame propagation, and thus continues to re- 
ceive attention for the insight it provides into these im- 
portant problems The problem was studied first by 
Saffman and Taylor (1958) 0, who injected air into oil in 
a Hele-Shaw cell and observed the formation of a single, 
steadily moving finger whose width decreased monotoni- 
cally to 1/2 of the channel width as the finger speed was 
increased. In subsequent experimental j^, numerical [Q, 
and theoretical work, the ratio of finger width to chan- 
nel width. A, was found to depend on a modified capillary 
number, 1/ B — 12 {w/hY Ca, which combines the aspect 
ratio, w/b, and the capillary number, Ca = (jV/c, where 
/i is the dynamic viscosity of the liquid, V is the veloc- 
ity of the tip of the finger, and a is the surface tension. 
For large l/B values, a transition to complex patterns of 
tip-splitting occurs m|,01. 

Our experiments reveal fluctuations in the width of the 
evolving viscous fingers that have not been reported in 
previous experiments |^,^,^ nor predicted theoretically. 
The fluctuations intermittently narrow the fingers from 
their expected width and are largest at low Ca, falling 
off as a power law with increasing Ca. Further, at large 
aspect ratios, the fluctuations are accompanied by sub- 
stantial departures from the monotonic dependence of 
the finger width on l/B found previously 1^-^: for as- 
pect ratios w/b > 250, which were not examined in pre- 
vious work, we find that the mean finger width no longer 
scales monotonically with l/B; for smaller aspect ratios, 
our finger width measurements are in accord with pre- 
vious results. Because large aspect ratios should more 
closely approach the ideal of two-dimensional flow, our 
observations pose a challenge to the assumptions under- 
lying theoretical analyses of viscous fingering. 



Experimental Methods. We conducted experiments in 
a 254 cm long channel formed of 1.9 cm thick glass plates. 
The spacers between the glass plates were stainless steel 
strips with thicknesses b = 0.051 cm, 0.064 cm, 0.102 cm 
or 0.127 cm; the channel width w between the spacers 
was varied between 19.9 cm and 25.1 cm. Both glass 
plates were supported at the sides so that they sagged 
similarly under their own weight (though by less than 
0.1% of the gap depth). For aspect ratios under 150, 
we used a smaller channel of length 102 cm and width 
7.4 cm. Interferometric measurements revealed that the 
root-mean-square variations in gap thickness were typi- 
cally 0.6% or less in the large channel and 0.8% in the 
small channel. Mechanical measurements of the bending 
of the glass due to the imposed pressure gradient revealed 
that such deflections were typically 0.2% or less; the max- 
imum deflection (2.2%) was measured in the widest chan- 
nel close to the oil reservoir at the highest flow rates. 
Experiments were conducted with air penetrating a Dow 
Corning silicone oil whose surface tension and dynamic 
viscosity was either a — 19.6 dyne/cm, fj, = 9.21 cP 
or (J — 20.6 dyne/cm, fj, — 50.8 cP at laboratory tem- 
perature (22°C). The oils wet the glass completely. A 
uniform flow rate was imposed by withdrawing oil with a 
syringe pump from a reservoir at one end of the channel; 
an air reservoir at atmospheric pressure was attached to 
the other end. 

The channel was illuminated from below and images 
were obtained using a camera and rotating mirror that 
captured 11 overlapping frames to produce concatenated 
images of up to 1200 x 10, 000 pixels at a resolution of 
0.25 mm/pixel. The interfaces were then digitally traced, 
yielding finger width values accurate to 0.1% in the larger 
channel and 0.3% in the smaller channel. For each flow 
rate, up to four time sequences of 20-30 digital inter- 
faces were recorded. Finger widths determined in con- 
secutive sequences agreed within the measurement ac- 
curacy. Mean width values agreed within 0.5% for data 
sets repeated after channel disassembly, cleaning, and re- 
assembly. 

For each interface image with a single, well developed 
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w/b = 158 w/b = 490 




1/B: 116 1757 408 2869 

V(cm/s): 0.0818 1.24 0.0301 0.211 

(Xy. 0.611 0.525 0.461 0.478 

5^: 1.30x10"^ 3.96x10"^ 3.32x10"^ 1.26x10"^ 

At(s): 270 18 810 108 



FIG. 1. Finger images recorded at regular time intervals 
for different values of aspect ratio w/b and modified capillary 
number 1/B, with corresponding values of the tip velocity V , 
the mean finger width (A), the rms fiuctuation of the finger 
width S\, and the time At between successive curves. 

finger, we found the instantaneous finger width by aver- 
aging measurements taken in a narrow window that was 
5% of the channel width beginning 1.2 channel widths 
behind the tip. This measurement window was chosen 
to be small compared to the length scale of the fluctu- 
ations (approximately a third of a channel width); the 
results were not sensitive to the window's exact location. 
From each time series of instantaneous widths, we de- 
termined the time average (A) and the root-mean-square 
(rms) fluctuation from the mean 6x. Each data set was 



analyzed for flow rates up to the point of tip splitting, 
beyond which the finger width A was no longer well de- 
fined II . 

Results. Typical interface image sequences are shown 
for w/b — 158 and 490 in Fig. |l|. For both aspect ratios 
the finger width A fluctuates visibly at low flow velocities 
(Fig. |l| (a),(c)). In the smaller aspect ratio system the 
width appears to become steady as the finger velocity 
is increased (Fig. |l|(b)), appearing exactly like the clas- 
sic "half-width finger" of Saffman and Taylor. However, 
with sufficient resolution, fluctuations can still be mea- 
sured for all velocities up to the onset of tip instabilities. 
In the higher aspect ratio system the width fluctuates 
visibly for all flow rates (Fig. ||(d)). The onset of tip in- 
stabilities in both cases occurs at 1/_B « 4000, similar to 
values seen in previous experiments |^,^ . 

For all the aspect ratios studied, we find that the rms 
fiuctuation of the finger widths is described by Sx = 
yl(Ca)'3 with A = (1.1±0.3) x IQ"" and /3 = -0.64±0.04, 
as Fig. H illustrates. We also observe that the instanta- 
neous velocity of the finger tip fluctuates from the aver- 
age velocity; within the experimental uncertainty these 
velocity fluctuations scale with Ca in the same manner 
as the width fluctuations. Interestingly, a dependence of 
the form Ga~'^^^ appears frequently in theories of viscous 
fingering 

The fluctuations in finger width are accompanied by a 
substantial deviation from the expected relation between 
finger width and velocity. Our results for the width of 
the viscous fingers for high w/b are not described by a 
single curve as a function of 1/B as predicted and 
the differences between data for different aspect ratios 
are far greater than those reported previously for low 
w/b Figure ||(a) illustrates this, showing the depen- 
dence of the mean finger width (A) on 1/B for values of 
w/b between 58.4 and 490. 

In particular, the mean finger width exhibits a surpris- 
ing maximum as the tip velocity is decreased for large 
aspect ratios (Fig. ^(a,b)). The value of (A) at the peak, 
('^)pcak' plotted vs. w/b in the inset of Fig. ||(a). 

To further compare our work with past results, we 
must account for the effects of the film of oil left be- 
hind on the plates as the finger advances. To do 
so, we use a correction to the modified capillary num- 
ber introduced by Tabeling and Libchaber, 1/B* — 
(1/S)/ [7r/4+ 1.7A(u;/6)(Ca)2/3] g. Using this correc- 
tion, they obtained agreement between their experimen- 
tal measurements of A and the theoretical predictions of 
McLean and Saffman g for 1/B* < 100 {1/B < 250); we 
obtain similar agreement for w/b = 58.4, as the inset of 
Fig. ||(b) illustrates. Fig. ||(b) shows the dependence of 
the mean finger width (A) on 1/B* for all aspect ratios. 
At high 1/B*, the film wetting correction also collapses 
the data to a common curve for all w/b > 58.4, though 
a slight displacement downward of the large w/b data 
remains. 
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FIG. 2. The rms fluctuation of finger widtli 5x as a function 
of capillary number Ca, where the dashed lines describe the 
best fit to all of the data sets within the observed scaling 
region, 5x = (1.1 x 10-*)Ca-° '^*. (The data for 394N and 
394W have the same aspect ratio but different widths; see 
caption of fig. 3.) The horizontal dash-dotted hnes in the top 
two graphs correspond to the limits of measurement accuracy 
for that channel and geometry; this limit is below the lower 
edge of the graphs for the other data sets. Fluctuations with 
5x ^ lO"'^ are not obvious visually. An upturn in 5x occurs 
at high Ca, signalling the onset of the secondary instabilities 
in the tips of the fingers. 

While our (A) results do not exhibit the classical scal- 
ing with 1/B*, the maximum value of the fluctuating fin- 
ger width observed during finger evolution, Amax, does, 
as shown in Fig. §(c). (Amax is the largest width value 
observed over time in a region behind the tip; observing 
only near the tip excludes widening due to relaxation ef- 
fects.) Amax does not exhibit a peak with decreasing 1/B, 
and while Amax has more statistical noise than (A), the 
Amax data collapse onto similar, monotonically decreasing 
curves that agree with the McLean-Saffman prediction at 
low 1/B* . The Amax data fall below the McLean-Saffman 
curve at high 1/B*, but not as much as the data for (A). 




200 400 600 

1/B* 

FIG. 3. (a) Mean finger width (A) vs. 1/B for various as- 
pect ratios w/b. Higher aspect ratios give smaller (A) values. 
The large w/b data exhibit a peak width ()^)pg^^^ as a function 
of l/B; the inset shows dependence of (A)pg_^]j on aspect ratio 
w/b. (b) Mean finger width (A) vs. 1/B* , which includes film 
wetting corrections. The solid line is the theoretical curve of 
McLean and Saffman [5]. The inset extends the w/b = 58.4 
curve to lower forcing, (c) Maximum width Amax of the fluc- 
tuating finger width observed during evolution as a function 
of 1/B* . Data for w/b > 137 collapse roughly onto the same 
curve, which follows the McLean-Saffman theory more closely 
than the low aspect ratio results. Aspect ratio symbols: o 58.4 
{w = 7.4 cm), X 137 {w = 7.0 cm), t> 158 {w = 20.1 cm), * 214 
(w = 22.6 cm), O 247 (w = 25.1 cm), + 314 {w = 19.9 cm), 
□ 394N {w = 20.0 cm), ★ 394W (it; = 25.0 cm), A 490 (w = 
24.9 cm). 
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These Amax data suggest that the fluctuations repre- 
sent an intermittent narrowing of the fingers from their 
"ideal" width. 

The finger width fluctuations and the peak in (A) ver- 
sus 1/B have proven robust under variations of experi- 
mental conditions. Both high and low viscosity oils gave 
the same results for the same geometric configuration. By 
treating the stainless steel spacers with an anti-wetting 
agent, we changed the contact angle where the interface 
is pinned at the back of the channel; this difference can be 
seen between Fig. ^(b) and(c); we again found the same 
results for the same geometric configuration. To ensure 
that the syringe pump was not imposing the observed 
velocity fluctuations, we pumped one data set by grav- 
ity siphoning, again with identical results. The effect of 
variations in the gap between the plates was examined for 
several aspect ratios by over-clamping the channel along 
the sides, which increases the intrinsic gap variations by 
a factor of 2.5 (measured interferometrically); both the 
finger width fiuctuations and the location of (A)pgjj)j. with 
respect to 1/B* were unchanged, although (A) values de- 
creased by about 4% near the peak. 

Though the fluctuation power law we observed re- 
mained unchanged for all experimental variations, we 
did discover that measurements on two channels with 
w/h = 394 but different values of w and h yielded dif- 
ferent results. The (A) values for both set-ups departed 
the common curve at approximately the same point with 
decreasing 1/B*, but at lower \/B* the behavior was 
different - compare 394W (394-widc) with 394N (394- 
narrow) in Fig. H(b). This difference suggests that either 
the width for a given w/b and Ca is not unique, or per- 
haps a third parameter, as yet unknown, is necessary to 
describe the problem for high w/h and/or low Ca. Con- 
sistent with this are the (A) data for w/h = 247 and 
314, which superpose closely even though they have no 
geometric parameters in common. 

It is unlikely that the fluctuations are caused by an in- 
stability of the film wetting layer because the film is very 
thin at low capillary numbers, where the fluctuations are 
largest. Film wetting fluctuations would also cause de- 
viations in the growth rate of the area of the fingers, 
which we do not observe. We speculate that the fiuctua- 
tions in finger width may be a consequence of long-time 
relaxations of the interface at the back of the channel, 
observable particularly in Fig. |l](c). We also speculate 
that the peak in (A) observed for decreasing 1/_B at large 
aspect ratios may also occur at small aspect ratios, but 
at values of too small to be reached. 

In conclusion, we have discovered fluctuations that in- 
termittently narrow evolving single fingers; the magni- 
tude of these fluctuations follows a power law with the 
capillary number for all aspect ratios studied. We also 
have found a departure from the classic scaling of flnger 
width versus 1/B for large aspect ratios [w/h > 250); 
the average finger width narrows at low 1/B, while the 



maximum finger width increases. These phenomena are 
not predicted by existing viscous fingering theories even 
though the phenomena are most pronounced for parame- 
ters that more closely match the theoretical assumptions. 
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